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ABSTRACT 

We have used the Green Bank Telescope to observe the millisecond pulsar 
O PSR J0514— 4002A on 43 occasions spread over 2 years. This 5-ms pulsar is 

located in the globular cluster NGC 1851; it belongs to a binary system and 
has a highly eccentric (e = 0.888) orbit. We have obtained a phase-coherent 
^ ; timing solution for this object, including very precise position, spin and orbital 

CS| ; parameters. The pulsar is located 4'/ 6 (about 1.3 core radii) from the center of the 

cluster, and is likely to lie on its more distant half. The non-detection of eclipses 
> ; at superior conjunction can be used, given the peculiar geometry of this system, to 

■ rule out the possibility of an extended companion. We have measured the rate of 

advance of periastron for this binary system to be Co = 0.01289(4)°yr _1 , which if 
due completely to general relativity, implies a total system mass of 2.453(14) M . 
Given the known mass function, the pulsar mass has to be <1.5 M and the mass 
of the companion has to be >0.96 M , implying that it is a heavy white dwarf. 
Oh" The 350-MHz flux density of this pulsar varies between 0.2 and 1.4 mJy; the 

O ' origin of these variations is not known. 

Subject headings: binaries: general — pulsars: general — pulsars: individual 
(PSR J0514-4002A) — globular clusters : general — globular clusters: individual 
'■ (NGC 1851) 

1. Introduction 

PSR J0514— 4002A was discovered in a 327-MHz search for steep-spectrum pulsars in 
globular clusters (GCs) carried out with the Giant Metrewave Radio Telescope (GMRT) at 
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Khodad, near Pune, India (Freire et al. 2004, henceforth Paper I). The pulsar is located in 
the globular cluster NGC 1851 and we will hereafter designate it to be NGC 1851A. It has 
a spin frequency of 200 Hz and is a member of a binary system, with a massive ~ 1 M 
companion and a very eccentric orbit (e = 0.888). 

In the Galaxy, all MSPs are either found isolated or in low-eccentricity binary systems 
with low- mass white dwarfs (WDs). There are pulsars with eccentric orbits and massive 
companions, but their spin periods are tens, not hundreds, of Hz. For these systems, there 
is not much time for accretion before the massive companion becomes a neutron star (see 
Lorimer, 2005 and references therein). 

NGC 1851A has hybrid characteristics. Its fast rotation indicates that the pulsar itself 
probably formed in the same manner as a normal MSP, ending up with a light WD companion 
in a nearly circular orbit. However, a subsequent stellar interaction could have disrupted 
the original binary system, replacing the lowest mass component (the original WD) with 
the present companion star in a highly eccentric orbit. Such exchange interactions are only 
likely to happen in environments with very high stellar densities, like the central regions 
of globular clusters. The same process could also form even more exotic systems, such as 
MSP/MSP or MSP/Black Hole binaries. 

An alternative scenario for the formation of such binaries is the nearly head-on collision 



of a neutron star (po ssibly even an MSP) with a giant star (e.g. iRasio &: Shapiro 1991 



Lombardi et al. 20061 ); this can form an e = 0.9 binary system consisting of the MSP and 
the stripped core of the giant star. We discuss this possibility at length in section §[5j since 
it might provide an explanation for the unusual flux density variations observed for this 
system. 

In Paper I, we had limited information to constrain the nature of the companion star, no 
mass estimates for either the pulsar or the companion, and no precise pulsar timing solution. 
In order to correct these deficiencies, we observed the pulsar on 43 occasions over the last 
two years with the Green Bank Telescope, and we present the results below. 



2. Observations and data reduction 
2.1. Data taking 

When observing NGC 1851A, we used mostly the prime-focus receiver PF1 with the 350- 
MHz feed. On several occasions, however, we used the PF1 receiver with the 425- and 820- 
MHz feeds, and twice used S-band receiver. The back-end instrument was the Pulsar Spigot 
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(Kaplan et al. 2005). For the observations at frequencies below 1 GHz, where the pulsar is 
clearly detectable, we observed a total bandwidth of 50 MHz, divided in 1024 channels, with 
a sampling time of 81.92 /xs. The S-band observations had a bandwidth of 600 MHz, 768 
channels, and 81.92 /is sampling. The data were dedispersed and folded using the PRESTO 
software packag^], resulting in 64-bin pulse profiles (~77/is per bin) approximately every 
500 s. 

With the pulsar's dispersion measure (DM) of 52.14pccm -3 , the 350 MHz observations 
had 492 [is of dispersive smearing in each channel, giving an effective system time resolution 
of 499 [is. At this frequency, dispersive smearing is the dominant contribution to the width 
of the Gaussian-like pulse profile, which is ~590 [is at half-power. At 425 MHz the dispersive 
smearing is 275 [is and the total time resolution is 287 [is, compared to the observed pulse 
width of ~370 [is. Therefore dispersive smearing is still the dominant contribution to the 
pulse width at 425 MHz. These measurements suggest that the use of coherent dedispersion 
would significantly improve the timing precision of this pulsar. 



2.2. Timing 



We cross-correla ted the average pulse profiles with a Gaussian template in the Fourier 
domain (ITaylor 19921 ) to obtain topocentric times of arrival ( TO As). These w ere then ana- 
lyzed with TEMPC^I, using the DE 405 Solar System ephemeris ( Standish 1998 ) to model the 
motion of the Green Bank Telescope relative to the Solar System Barycenter. 

After an initial fit, we obtained a x 2 P er degree of freedom of 1.5. Because there are 
no unmodeled trends in the residuals (their distribution is Gaussian), we believe that the 
routine used to derive the TOAs is slightly under-estimating the TOA uncertainties. After 
increasing the TOA uncertainties by 22% (a value similar to what we have seen for isolated 
MSPs timed with the same software) we obtained a \ 2 per degree of freedom of 1.0. The 
resulting timing parameters are presented in Table [21 the uncertainties are the 1-a TEMPO 
estimates. The significance of these timing parameters is discussed below. 

The orbital paramete r s were determined using the Damour & Deruelle orbital model 
( IDamour fc Deruelle 19851 ; iDamour fc Deruelle 19861 ) . We measured the longitude of peri- 
astron (u) and the time of periastron passage (T ) at a time (August 2005) where we have a 
sequence of daily observations that cover one full orbit. That T was chosen as the reference 
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epoch. The relativistic parameter 7, if included in the model, is very strongly correlated with 
the projected semi-major axis x and uj. When 7 eventually becomes needed by the timing 
model, the measured x and uj will not change by more than 0.007s and 0.09°, respectively, 
from their present values. 



3. Position, period and period derivative 

In Paper I, the position of the source identified as the pulsar in an interferometric map 
of the cluster is given as a = 05 h 14 m 06.74±0.06 s , S = -40°02'50.0±1.3" (J2000). The a and 
5 given in Table [2] are l-a consistent with these values, but much more precise. The pulsar 
is about 4'.'6 (about 1.3 core radii) East of the center of the cluster (see Tabled]). Such close 
proximity to the center is typical for GC pulsarg§ due to the effects of mass segregation. The 
uncertainty in the derived magnitude of the projected offset from the center of the cluster is 
dominated by uncertainties in the position of the cluster center and the cluster's core radius. 



The observed period derivative, P bs, is the sum of several different terms (jPhinney 19921 ): 




obs 




- (a G + ii 2 D + a c ) 

c v ' 



where P; nt is the intrinsic period derivative of the pulsar, a G is the difference in Galactic 
accelerations between the Solar System Barycenter and NGC 1851 projected along the line- 
of-sight, n and D are the proper motion and distance to the pulsar (assumed to be the 
same as the cluster values in Table [T]) and a c is the acceleration of the binary caused by 
the cluster, also projecte d along the line-of-sight. For a G , we obtain 6.1 x 10~ 12 ms~ 2 



( IKuijken fc Gilmore 19891 ) . For fi D, we obtain 6.5 x 10 ms , using the value for fx 



obtained by Dinescu et al. (1999). With these values and P D b s we obtain: 

P c = P int + P^ = (+0.0 ± 1.4) x 10- 22 . (2) 

c 

Pulsars with negative P c can put lower limits on the absolute value of a c , and for that 
reason have been used to study the gravitational potential of their parent GCs (e.g., Possenti 
et al. 2002, Freire et al. 2003). In this case, that can't be done, but we can say that it is 
very likely that the pulsar is located on the more distant half of the cluster. If it were on the 
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near side, then a c in eqn. [2] would be positive, and Pi nt < 2.8 x 1CT 22 (2— o upper limit). 
This would imply an anomalous characteristic age of more than 280 Gyr. In the likely event 
that the pulsar has a more normal characteristic age, then a c has to be negative, placing the 
pulsar on the far side of the cluster. 

In Tabled! v z (0) is the spread of stellar line-of-sight velocities at the core. Using this, 
the cluster core ra dius and an analytical gravitational model for the central regions of GCs 



( jFreire et al. 20051 ). we can calculate the maximum negative acceleration that the cluster 
can induce on the pulsar at its sky position. From this we can derive an upper limit on 
the intrinsic period derivative and surface magnetic flux density, and a lower limit on the 
characteristic age (see Table [2]), as assumed in the formation scenarios in Paper I. 



4. On the nature of the companion 

4.1. Is the companion extended? 

In Paper I, we reached the conclusion that, in the event of an extended companion, the 
circularization timescale for this system is larger than the age of the cluster. An extended 
companion could not therefore be disproved because of the system's eccentricity. 

In Table |3l we present results of calculations of the companion masses, absolute com- 
ponent separation at superior conjunction (a sup ) and its projection in the plane of the sky 
(S = a sup cosz) for a set of orbital inclinations. The companion masses are calculated for a 
pulsar mass of 1.2 M ; the mass of the lightest known NS, the companion of PSR J1756— 2251 
(Faulkner et al. 2005, Ferdman et al. 2006). For these companion masses, we list the radii of 
zero-age main-sequence stars, calculated using log in (-R/- R^) = 0.917 x \o^ n ( M/M m ) — 0.020, 



which apply for stars of mass similar to that of the Sun (lAllen fc Cox 20001 ) . The use of the 
smallest possible pulsar mass and corresponding zero age main sequence radii leads to the 
lowest possible companion masses and (if extended) sizes, the most conservative assumption 
in the discussion that follows. 

For i < 60° (limits are approximate), the companion masses are above 1M Q . In a 
star cluster as old as NGC 1851 (about 9 Gyr, see Tab. [TJ, such stars have either finished 
their life cycles and left compact remnants behind, or are now in their giant phases. Giant 
stars have radii that are tens to hundreds of times that of similarly massive stars on the 
main sequence, which are of the order of one solar radius. Such giant stars would not even 
fit within the distance between the components of NGC 1851A at superior conjunction or 
periastron. In the first scan taken in 2006 July 7, we can detect the pulsar clearly through 
periastron, implying that the possibility of such a giant companion is excluded. 
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For i > 80°, the radius of a main- sequence companion star is larger than S, implying 
that an eclipse should occur. The longitude of periastron for this system is 82.27°, and 
therefore superior conjunction occurs at a true anomaly of 7.73°, 15min and 54.1s after 
periastron. For the second scan taken on 2006 July 7, the two TOAs occur 11.3 and 22.5 
min of barycentric time after the periastron^. The pulsar is detectable throughout the whole 
scan. This excludes the possibility of the companion being a main sequence star with a mass 
smaller than 0.87 M Q . 

For 60° < i < 80°, we could in principle avoid eclipses and still postulate a main- 
sequence companion. For any higher pulsar mass, or extended stars larger than zero age 
main-sequence, this inclination interval becomes smaller. However, even for the lowest incli- 
nation listed in Table [3], the pulsar signal would be passing only 1.2 R Q above the surface of 
the companion, right through its corona. We also note that the pulsar wind should lead to a 
denser and more extended corona in this star, so we should be able to detect some increase 
in the electron column density at this orbital phase. Dividing the second observation of 2006 
July 7 in 4 sub-bands and calculating TOAs for the whole observation, we obtain a DM of 
52.148(5) pccm -3 , not measurably different than the DM at other orbital phases and smaller 
than the increases in the plasma column density seen at superior conjunction for pulsars with 
extended companions (see Freire, 2005 and references therein). 

Furthermore, if the companion was extended, there would be tidal effects on the orbit 
of the pulsar (e.g., PSR J0045-7319, Lai 1996 and PSR B1259-63, Wang et al. 2004). 
These include changes in the inclination (and therefore, changes in the projected size x = 
a p /cd(smi)/dt) of the orbit, large variations in the orbital period (Pb ^ 0), and unmodeled 
and systematic trends in the residuals. As described above, the values for x and Pb are not 
significant (see Table [2]), and the object times as well as an isolated MSP, with no unmodeled 
trends in the residuals despite the fact that our timing precision is much higher than for either 
PSR J0045— 7319 or PSR B1259— 63. In order to avoid eclipses and any measurable tidal 
effects, the companion has to be a compact object. The implications of this are discussed 
below. 



4 Because of the advance of periastron between the epoch of the parameters in Tableland this observation, 
superior conjunction occurred 15 min and 52.8 s after periastron, a difference of 1.3 s relative to the reference 
epoch. 
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4.2. What is the mass of the companion? 

For a system consisting of two compact objects, the observed u is due solely to the effects 
of general relativity. This allows an estimate of the total mass of the system: 2.453(14) M . 
Given the mass function, the pulsar mass cannot be larger than 1.50 M and the companion 
mass must be larger than 0.96 M . For a median inclination of 60°, the mass of the pulsar 
is 1.350 M Q , a value that is fairly typical of the neutron stars with well-determined masses. 
In this case, the companion mass would be about 1.105 M Q . 

For 49.76° < i < 52.86°, both components would have masses within the present range 
of well- measured neutron stars (see Fig. Q]); however, given a flat probability distribution in 
cosi, it is about 15 times more likely that 52.86° < i < 90°, where m c < 1.2 M . In the 
higher inclination range, there is the possibility that the companion is the lightest NS ever 
discovered. However, this is probably more than compensated for in the lower inclination 
range by the possibility that the companion is a WD or a stripped core of a giant star with 
m c > 1.2 M . We therefore believe that the probabilities are not far from 15 to 1 against the 
companion being a neutron star. Given the slight possibility, though, that the companion 
star could be a MSP, we searched several of the 350, 820, and 1950 MHz observations for 
additional pulsations, but found none. 

5. Flux Density Variations 

In Table|5]we list the average flux densities at 350, 425, 820 and 1950 MHz, together with 
the number of observations averaged, we derive a spectral index of —2.65. At 350 MHz, the 
flux density changes very significantly, from about 0.2 to 1.4 mJy. Even our brightest detec- 
tion is significantly fainter than the average flux density reported in Paper I (3.4 ± 0.4 mJy). 
This could explain the discrepancy in reported spectral indices; in Paper I we use the non- 
detection at 610 MHz to derive as < —3.4. However, the GMRT flux density at 325 MHz 
comes from imaging, so it is possible that there is extra unpulsed radio emission in the 
vicinity of the pulsar. However, we can not exclude the possibility of some lapse in our 
understanding of the sensitivity of these systems. 

The timescale for the variations near periastron (r p ) is of the order of 1000 s; near 
apastron the flux density is seen to change very slowly and monotonically over several hours, 
only becoming significant from one day to the next (i.e. T a /r p > 10). 

The flux density variations are almost certainly not caused by diffractive scintillation. 
The scintillation bandwidth at 350 MHz is (very roughly) 10— 20 kHz (Cordes and Lazio 
2002), so our observations average over many hundreds of scintles in the 50 MHz band- 
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width. Furthermore, refractive scintillation will change the flux by no more than ~50%, on 
timescales of ~10 7 seconds (Cordes and Lazio, 2002). Therefore, if the ISM between us and 
the cluster has everywhere a Kolmogorov spatial spectrum and is described, to within a few 
orders of magnitude, by the Cordes and Lazio model, then no large-amplitude variations 
in the flux density should be observed. Scintillation by a region or regions of particularly 
high ionized gas density could provide an explanation, we could alternatively be dealing with 
small eclipses caused by gas in the vicinity of the binary. 



5.1. Gas from the Companion Star 

Flux density variations of similar magnitude and duration have been observed for several 
eclipsing binary pulsars away from superior conjunction. The best documented examples are 
PSR Jl 740— 5340, a MSP-main sequence binary located in the GC NGC 6397 (Ferraro et 
al. 2001), and PSR J1748-2446A , an eclipsing binary MSP located in the GC Terzan 5 



( ILyne et al. 19901 ; iLyne et al. 20001 ) also known as Terzan 5 A. As in virtually all eclipsing 
binaries, the companions of these two pulsars are losing g result of their interaction 

with the pulsar wind. Part of the observed flux density variations away from superior 
conjunction are due to secondary eclipses thought to be caused by gas clumps that remain in 
the vicinity of these binary systems; even though the precise mechanism of the modulation 
is not known. Many, but not all, of these secondary eclipses are "total" (i.e. the pulsars 
become completely undetectable). 

The companion of NGC 1851A is much more massive than the companions of either 
PSR J1740— 5340 or Terzan 5 A. No outgassing has ever been clearly detected by pulsar 
timing in systems where the companion is a massive white dwarf, and it has not been 
detected in NGC 1851 A either (Q. However, there could still be gas clumps in the vicinity 
of the system if the companion of NGC 1851A is outgassing at a lower level than what we 
can detect. This is a possibility if it is a stripped core of a giant star. 

Unlike PSR J1740-5340 or Terzan 5 A, in NGC 1851A we detect no "total" eclipses 
caused by gas clumps at superior conjunction or away from it. Furthermore, a common 
feature of such eclipses is that they invariably affect the lower frequencies more strongly 
(e.g. longer duration eclipses), although in varying degrees for different binary systems. 
That is not observed in NGC 1851A, where the attenuation may be frequency dependent, 
although in a variable and inconsistent manner. However, it is important to note that the 
behavior of the other pulsars mentioned above varies strongly from orbit to orbit. In our 
case, we only have good observations through a single periastron and superior conjunction. 
With the current data we cannot conclusively demonstrate that the variations are caused by 
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orbiting gas. However, we do consider it a possible explanation for the observed flux density 
variations. 



5.2. Refractive scintillation 

It is possible that, despite the predictions of Cordes and Lazio, the changes in flux 
density are indeed due to scintillation. This is more likely if, somewhere between us and the 
pulsar, there are unusual regions of the ISM acting as strong lenses. We can use the changes 
in t to constrain their location. Assuming a maximum inclination of 90°, the pulsar is 
traveling at v p ^lTSkms" 1 relative to the center of mass at periastron, and at v a ~10kms _1 
relative to the center of mass at apastron. Furthermore, the distance and proper motion 
listed for this cluster in Table [1] imply a perpendicular cluster velocity of ~160 km s _1 . 

If the lenses are near us, at apastron the velocity of the pulsar v' a relative to those 
lenses can not be much different from the perpendicular velocity of the cluster times / (the 
distance from the Solar System to the lens divided by the distance from the Solar System 
to NGC 1851). At periastron, we obtain (20 x / < v' p < 320 x /) km s -1 , depending on 
how the pulsar's orbital velocity at periastron adds to its proper motion. This would imply 
0.1 < v' p /v' a = T a /r p < 2, which is inconsistent with observation. If / is very small, then 
the Earth's orbital velocity becomes dominant in determining r and we can then produce a 
larger range of T a /r p , however, no seasonal variations in r a are observable. If, on the other 
hand, the medium responsible for the variations has a small velocity relative to the center of 
mass of the binary, we obtain r a / r p ~ fv p / fv a = (1 + e)/(l — e) = 16.85, which is consistent 
with observation. 

This suggests that a hypothetical lensing medium would probably be located in the 
cluster, which is also the case if we were instead observing gas clumps orbiting the binary 
itself. Given the proximity to the pulsar, we can then readily estimate the associated spatial 
scale of the medium, independently of its nature, from v p and r p . about 10 5 km. 

6. Conclusion and prospects 

In this work we have determined the phase-coherent timing solution of NGC 1851A. This 
allowed us to locate the pulsar in the globular cluster and determine precise rotational and 
orbital parameters. The geometry of the system, the lack of any observable tidal effects in the 
timing and the lack of eclipses or a measurable increase of the DM at superior conjunction 
allow us to exclude the possibility of an extended companion. This in turn allows us to 
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calculate the total mass of the system from the observed rate of the advance of periastron: 
2.453(14) M . The pulsar mass has to be smaller than 1.5 M and the companion mass has 
to be larger than 0.96 M . The companion is either a massive WD or the stripped core of a 
giant star, and less likely a neutron star. 

The flux density variations remain puzzling. They could be due to gas clumps in the 
vicinity of the binary created by outgassing from the companion, a possibility considering 
that it could be the stripped core of a giant star. Alternatively, they could be due to refractive 
scintillation, in which case the lensing structures are also likely to be located in NGC 1851. 
his would imply the presence of intra-cluster gas, the second instance after 47 Tucanae 



(IFreire et al. 20011 ) where such gas was confirmed in a globular cluster. 



If we observe the pulsar often close to periastron, it might be possible to measure 
the velocity of the binary relative to the medium that causes the flux density variations 
and decide which of the scenarios discussed is the correct one. Furthermore, a detailed 
study of the orbital variability of r might allow a determination of the orbital inclination 
of the system, and lead to a determination of the individu al masses. So far, o nly one MSP, 



PSR J1909-3744 has a well-determined mass, 1.44 M ( IJacobv et al. 20051 ) . Measuring 
more MSP masses precisely is important in order to estimate how much mass a neutron star 
needs to accrete in order to become a MSP. 

Such a detailed mass measurement will be confirmed (or refuted) by the measurement 
of another relativistic timing effect, the Einstein delay (7). Our simulations suggest that, 
with the present observing strategy and timing precision, we will be able to measure 7 to 
better than 1 ms within 6 to 8 years. For an edge-on orbit, 7 is predicted to be 16.13 ms, and 
for % = 60°, 7 = 19.44ms. The Shapiro delay can not be measured with the present timing 
precision, no tests of general relativity are possible with the present telescope sensitivities. 
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paper. 
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Table 1. Parameters for the globular cluster NGC 1851 



Right Ascension of center, a 


05 h 14 m 06!3 


(Harris 1996) 


Declination of center, 5 


-40°02'50 


(Harris 1996) 


Galactic longitude, I (°) 


244.51 




Galactic latitude, b (°) 


-35.04 




Cluster distance, D (kpc) 


12.1 


(Harris 1996) 


Core radius, 8 C (arcmin) 


0.06 


(Harris 1996) 


v z {0) (km s" 1 ) 


11.3 


(Dubath, Meyland & Mayor 1997) 


Proper Motion in a, \x a (mas yr -1 ) 


+ 1.28 ±0.68 


(Dinescu et al. 1999) 


Proper Motion in 5, fig (mas yr _1 ) 


+2.39 ±0.65 


(Dinescu et al. 1999) 


Age (Gyr) 


8.9±0.9 


(Salaris & Weiss 1997) 
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Table 2. Parameters for the PSR J0514-4002A binary system 



Observation parameters and flux density 



First observation (MJD) 


53258 


Last observation (MJD) 


54066 


Number of TOAs 


316 


Residual rms (/zs) 


30 


Epoch (MJD) 


53623.155088 


Flux density, 




... at 350 MHz, S350 ( m Jy) 


0.62(2) (93) a 


... at 425 MHz, S 425 (mJy) 


0.28(6) (2) 


... at 820 MHz, S S 20 (mJy) 


0.065(6) (12) 


... at 1950 MHz, S1950 (mJy) 


0.0056(10) (3) 


Spectral index, as 


-2.71(9) 


Timing parameters 


Right Ascension, a (J2000) 


05 h 14 m 06?6927(2) 


Declination, S (J2000) 


-40°02'48'.'897(2) 


Spin frequency, v (Hz) 


200.37770740529(11) 


Time derivative of v, v (10~ 17 Hz 2 ) 


-4.7(5) 


Dispersion Measure, DM (pccm~ 3 ) 


52.1489(6) 


Orbital period, P\, (days) 


18.78517915(4) 


Projected size or orbit, x (1-s) 


36.296588(9) 


Time of passage through periastron, T (MJD) 


53623.1550879(4) 


Eccentricity, e 


0.8879773(3) 


Longitude of periastron, ui (°) 


82.266550(18) 


Rate of advance of periastron, u (° yr _1 ) 


0.01289(4) 


Einstein delay, 7 (ms) 


[8±26] b 


Time derivative of x, x (10 -12 1-s/s) 


[-0.1±0.4] 


Time derivative of P B , Pb (10^ 9 ) 


[1.5±0.7] 


Derived parameters 


Angular distance from center of cluster, 6± 




... (arcseconds) 


4.6 


... (core radii) 


1.3 


... (pc) 


0.27 


Spin period, P (ms) 


4.990575114114(3) 


Time derivative of P, P (10~ 21 ) 


1.17(14) 


Maximum cluster acceleration, a Cima x (m.s -2 ) 


18.4 x 10~ 9 


Maximum intrinsic P, Ant, max (10 -19 ) 


3.0 


Maximum magnetic flux density, Bq (10 9 Gauss) 


1.2 


Minimum characteristic age, r Cim i n (Gyr) 


0.26 


Mass function, / (M Q ) 


0.14549547(11) 


Total system mass, M (M Q ) 


2.453(14) 
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Table 2 — Continued 

Maximum pulsar mass, Mj, imax (M©) 1.50 
Minimum companion mass, M Ciinin (M Q ) 0.96 

a Number of measurements averaged in parenthe- 
ses. Flux uncertainty is 30% / VN 

b Values in square brackets arc not 3-cr signifi- 
cant, nor are they expected to be so. They were 
not fit when determining the remaining timing pa- 
rameters. 
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Table 3. Geometric parameters for the companion of NGC 1851A 



i 


m c 


&sup 


S = a S up cos i 


R c 


(°) 


(M ) 


(Ro) 


(R Q ) 


(Ro) 


90 


0.8483 


4.2471 





0.82 


85 


0.8528 


4.2502 


0.3704 


0.83 


80 


0.8665 


4.2596 


0.7397 


0.84 


75 


0.8901 


4.2758 


1.1067 


0.86 


70 


0.9252 


4.2996 


1.4705 


0.89 


65 


0.9739 


4.3322 


1.8309 


0.93 


60 


1.0396 


4.3754 


2.1877 


0.99 




Fig. 1. — Constraints on the masses of NGC 1851 A and its companion. The hatched region 
is excluded by knowledge of the mass function and by sin % < 1 . The diagonal dashed lines 
correspond to a total system mass that causes a general-relativistic Co equal or within 2-er of 
the measured value. The four solid curves indicate constant inclinations. We also display the 
probability density function for the mass of the pulsar (top) and the mass of the companion 
(right), and mark the respective medians with vertical (horizontal) lines. 



